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Abstract 

We have measured the low frequency critical current noise in Nb/AlOx/Nb Josephson junctions. 
Unshunted junctions biased above the gap voltage and resistively shunted junctions biased near 
the critical current, /c, have been measured. For both, the spectral density of SIc/Ic, Si^{f), is 
proportional to 1//, scales inversely as the area. A, and is independent of Jc = Ic/A over a factor 
of nearly 20 in J^. For all devices measured at 4.2 K, 5*4^(1 Hz)= 2.0 it 0.4 • W~^^ /Hz when scaled 
to yl = 1 /im^. We find that, from 4.2 K to 0.46 K, Si^{f) decreases linearly with temperature. 
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There has been a renewed interest in understanding the ubiquitous 1// fluctuations in 
the critical current and resistance in Josephson junctions. This is motivated in part by the 
intense effort to reduce the noise-induced decoherence in a variety of qubits based on these 
junctions. The dominant model for these fluctuations is a collection of two level fluctuators 
(TLFs) in the oxide barriers of the junction that randomly and independently switch on and 
off the conductivity and supercurrent in small localized areas of the junction. This model 
P] gives rise to the expectation that the normalized spectral density of the critical current 
and resistance noise, e.g. Si^ = Sjc/I^ or Srj = Srj/Rj, should be roughly independent of 
the junction conductivity (or critical current density, Jc) and should scale inversely as the 
junction area, A. It has been shown that the IcRj product of the junction is unaffected by 
these fluctuations Thus the two normalized spectra should be equal. 

A survey of critical current fluctuations in Nb junctions from a number of laboratories 
{2! recently found that this normalized spectral density, scaled to an area A = 1 /im^, was 
indeed roughly constant, had a value of 1.4- 10~^'^/Hz at 4.2 K and varied with temperatures 
as T^. These results were then explained theoretically More recently, measurements of 
resistance fluctuations [4] in submicron Al junctions made using shadow evaporation gave 
rather different results: the normalized and scaled spectral density was nearly three orders 
of magnitude smaller at 4.2 K, but varied linearly with temperature from about 100 K 
to nearly 0.8 K, where it appeared to saturate. It has been suggested that a possible 
explanation for this is that fluctuations in normal and superconducting junctions arise from 
different mechanisms. Very recently a prediction has appeared jsl , based on the independent 
fluctuator model, that these fluctuations should vary linearly with T but should increase as 
the flfth power of the thickness of the barrier, t^. Since both Rj and Ic depend exponentially 
on tb, this implies a small, but potentially detectable, dependence of the normalized spectral 
density of the fluctuations of these parameters on their values. In this Letter, we present 
data for the dependencies of these fluctuations on junction area, critical current density 
and temperature in Josephson junctions fabricated using a well-characterized, very versatile 
trilayer (Nb/AlOx/Nb) process where the values of the barrier thickness as a function of Jc 
have been accurately determined. 

Examples of the normalized noise spectra we observe from our junctions are shown in 
Fig. [H Here and throughout the paper, measured noise spectra are presented in the terms 
of an equivalent current noise source. Si (or Si = Sj/I^), shunting the junction. Spectra 
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shown are for three junctions with areas ranging from 3 fim^ to 94 yum^ measured at 4.2K. 
These junctions, having Jc = 100 A/cm^, have no external shunts and are biased at 3.7 
mV, which is above the gap voltage. The data were obtained by placing the junction as 
one arm of a Wheatstone bridge having a commercial SQUID as a null detector. In each 
case the spectra are fit to the sum of white noise plus a 1// spectrum, 5*1//. The white 
noise component is just that expected from the shot noise due to the bias current, which 
dominates other white noise sources, e.g. the Johnson noise due to the other arms of the 
bridge circuit. Si^ is obtained from the 1// component of the spectra corrected for the bias 
conditions, e.g. for unshunted junctions biased above the gap voltage the 1// noise should 
almost be entirely due to resistance fluctuations. For this case, the constancy of the IcRj 
product gives Si^ = Si/f^. 

Throughout the paper, the values presented for 5*1// and related quantities are at 1 Hz. 
The inset of Fig. [1] shows 5*1// for the three junctions as a function of A, along with a 
line having the predicted linear dependence with A. This clearly demonstrates the scaling 
expected from the independent fluctuator model. To facilitate the comparison of Si^ of a 
variety of junctions, we present data for Si^ scaled to A = 1 yum^ and denoted by S-^. For the 
junctions in Fig. [H S-^ = 2.2 ■ lO^^^/Hz, which is about two orders of magnitude less than 
that commonly seen in other superconducting junctions We have further verified that 
the noise is consistent with resistance fluctuations by observing that 5*1// is proportional 
to /| between 3.7 mV and 6 mV, i.e. Sjij is independent of lb, where lb is the total dc 
current through the junction. Such unshunted junctions are not suitable for a study of 
the temperature dependence of this noise due to the high power dissipation at the required 
bias voltage that leads to substantial electron and lattice heating at lower temperatures [3] ■ 
Therefore, to extend this work to lower temperatures we fabricated shunt resistors in parallel 
with the junction to form a non-hysteretic resistively shunted junction (RSJ). 

All junctions reported here were^abricated in our laboratory using a process (SAL-EBL) 
that has been detailed previously [8]. These junction have been shown to be of very high 
quality and dimensional uniformity. This high quality is also demonstrated by very small 
subgap leakage currents, e.g. a subgap resistance, measured at 400 mK and 0.5 mV that is 
over 10^i?j The initial step in the fabrication of the Nb/AlOx/Nb junctions is the in 
situ sequential deposition of the trilayer through a liftoff resist mask. The tunnel barrier is 
formed by thermal oxidation of the 80 nm Al interlayer in dry O2 for exposures of 10^ — 10^ 



3 



Torr-min (10^ — 10^ Pa-s) giving the range of Jc of 50 — 800 A/cm^ studied here. The 
properties of the barriers, i.e. barrier height and thickness, have been carefully determined 
by fitting their normal conductance vs. voltage using numerical solutions of the Schrodinger 
equation [9j. 

For RSJs, the shunt resistors (see Fig. [21 inset) were designed to minimize their parasitic 
inductance, e.g. by placing a Nb ground plane under the resistor, while having a sufficient 
volume and contact area to ensure adequate cooling of both the electrons and the lattice. 
To ensure nonhysteretic I-V curves, the resistance of the shunt resistor was selected to give 
(5c — 0.3, where /3c = ^iielcR^C /h. With this constraint the ratio of the 1// noise to 
the Johnson noise due to the resistor should scale as Jc^^. This fact, along with rapidly 
increasing oxidation times, sets a practical lower limit on J^- The upper limit is set by a 
de grad ation of the barrier quality, as seen by increasing subgap leakage, for Jc > 5 kA/cm^ 

The temperature rise, ST, due to the power, P, dissipated in the resistor can be described 
by 6T = {kP + To")i/" - To, where k is a constant depending on the material and geometry 
and To is the substrate temperature. We have measured this heating in resistors having 
the same design as our shunt resistors by observing the increase in resistor temperature 
(measured using noise thermometry) as a function of P over the range of temperatures 
used. These data are shown in Fig. [2] for temperatures ranging from 0.35 K through 0.7 
K. As one can see, the assumed fitting functions with n = 4.1 ± 0.1, which is consistent 
with lattice heating [3], and k = 22.5 ■ lO^'^ ± 0.9 (K"/nW) provide good descriptions of the 
temperature increase, and can be used to provide heating corrections to the data if needed. 
The noise spectra for the RSJs were generally measured at an input power of about 1 nW. 
For this low power, heating makes a significant, though small, correction to the temperature 
of the RSJs at the lowest temperatures in the measurement. 

The junctions were cooled using a ^He refrigerator having a base temperature of 350 mK. 
Filtering at room temperature and 4.2 K along with battery powered amplifiers and bias 
current sources were used to isolate the junctions from electronic noise. This filtering, along 
with the bridge measurement circuit reduced the bias current noise to levels too low to be 
detected. The upper arms of the bridge are ~ 1 kQ resistors anchored at 4.2 K. The lower 
arms consist of the sample along with a matched balancing resistor, Rbai- This balancing 
resistor was either a metal film resistor attached to the chip carrier in the case of unshunted 
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junctions, or, for RS Js, an on-chip resistor of the same design as the junction shunt resistor, 
Rsh (see Fig. 2 inset). The SQUID null detector had a noise limit of 0.5 pA^/Hz above 10 
Hz and a 1// knee at 1 Hz when directly connected to the bridge. 

The noise spectra from the RSJs could also be fit to the sum of white noise plus 1// 
noise. In general, the values of these two quantities were anomalous at higher voltages where 
clear resonances were observed. These behaviors can be seen in the inset of Fig. [3l which 
shows the dynamic resistance along with 5*^*^ and 5*1(1 kHz). These quantities are plotted vs 
Ai = (/ft — Jc) //c on a log scale so the behavior near 1^ is clearly visible. The resonances were 
due in part to the parasitic inductance {Lp = 2.3 pH) of the shunt resistor. This behavior 
made it impossible to use the value of the white noise as a reliable thermometer. However, 
the values of Si^ extracted from the spectra are independent of bias at low voltages, i.e. 
for bias currents just above Ic- In this bias-independent region, the values of Si^ obtained 
also agrees with those obtained from unshunted junctions on the same wafer at 4.2 K where 
heating in the unshunted junctions is not an issue. Figure [3] compares the values of 5*^*^ 
obtained at 4.2 K for both shunted and unshunted junctions with areas in the range of 3 
fim^ to 94 fim^ and Jc values from 55 to 807 A/cm^. Within the scatter of about a factor 
of 2, 5*^*^ is independent of area and Jc as expected from the simple independent fluctuator 
model. The temperature dependence of S-^ for our RSJs measured near Ic is shown in Fig. 
m A best fit (dashed line) having a linear temperature dependence clearly describes the 
data quite well. At 4.2 K SI = 1.42 • IO^^VHz and decreases to 0.16 • IO^^VHz at the lowest 
temperature, 0.46 K. The lowest temperature data points were corrected for heating in the 
shunt resistor at the measurement power of about 1 nW. The largest correction at the lowest 
temperature was ~90 mK, but the conclusion that S-^ varies linearly with T would not be 
significantly changed even if this correction were omitted. 

A recent analysis of critical current noise due to TLFs in the junction barrier ^ predicts 
that this noise is given by S^^ ~ (T / f)t\, i.e. it increases as the fifth power of the junction 
barrier thickness, tb- The solid line in Fig. [3] shows the best fit to this equation using the 
proportionality constant as the fitting parameter. This proportionality constant contains 
the density of states of the TLFs, which is therefore assumed to be independent of Jc- We 
have previously done careful comparisons of the dependencies of both Jc and ti, on oxidation 
parameters , which allows us to accurately determine the barrier thickness corresponding 
to a given Jc- While the predicted variation of Sf on is not seen, we must note that there 
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is no independent measure of the density of TLFs. There is some reason to suspect that this 
density is increasing with Jc, since for Jc > 5 kA/cm^ it is well known that barrier leakage 
increases significantly [lo| . 

In summary, we have measured 1// fluctuations in the critical current and resistance of 
Josephson junctions and observed their dependence on junction area, critical current density 
and temperature. We find that Sj^ scales linearly with temperature, without saturation, 
down to the lowest temperature measured (0.46 K). In addition, S^^, averaged over all 
junctions measured, gives 2.0 ± 0.4 • 10~^^/IIz at 1 Hz and 4.2 K, which is nearly two orders 
of magnitude less than that commonly observed in Josephson junctions. However, the data 
do not exhibit the recently predicted dependence of S-^ on barrier thickness as t^. 

This work was supported in part by NSF grant #DMR-0325551. 
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FIG. 1: Current noise spectra, Sj, of unshunted junctions having areas of 3 /xm^, 22 ^um^ and 94 
fim^ at 4.2 K and Jc = 100 A/cm^. The junctions were biased at 3.7 mV. The sohd curves through 
the spectra are best fits to white plus 1/f noise as described in the text. The inset shows the 1/f 
component of the noise spectra at 1 Hz as a function of junction area with a hnear fit to the data 
(red hne). 

FIG. 2: Temperature increase as a function of power dissipated in AuPd thin film resistors for bath 
temperatures of 0.35 K to 0.7 K. The length (L), width (W) and thickness (t) of each resistor are 
L = 33.5 fMui, W = 28 /xm, t = 0.08 /xm. The inset shows the Wheatstone bridge circuit with a 
thin film AuPd resistor and an RSJ in the lower half of the bridge circuit (note that the RSJ was 
replaced by a sample resistor for these data). 

FIG. 3: Sf^, at 4.2 K of shunted (solid black symbols) and unshunted (open blue symbols) junctions 
as a function of Jc for nominal junction areas of 2.25 ^m^ (triangles), 3 fim^ (circles), 22 fim^ 
(diamonds), 94 /xm^ (squares). The solid red line is a best fit to the data from ref. {o^ as described 
in the text, and the dashed line is the value from ref. [2]. The inset shows Sf^ (solid black circles) 
and the normalized noise at 1 kHz, 5*1(1 kHz), (open green triangles) versus Ai for an RSJ having 
Jc = 807 A/cm^ and A = 3.3 ^m . The green dashed line in the inset is the expected value of Sf(l 
kHz). 

FIG. 4: Temperature scaling of normalized critical current spectral density at 1 Hz for two 3.3 
fim^ RSJs from two different process runs represented by closed and open symbols. The RSJs 
had Ic = 27 /xA and R^h = 4.4 0,. The dashed line shows the best fit to the data having a linear 
dependence on T. 
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